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ABSTRACT 


All  the  radiative  corrections  are  calculated  to  second  order 
in  elementary  charge  e  for  the  decay  K—  jj~  +  rr°  +  V  with  finite  neu- 
tretto  mass.  The  infrared  divergent  terms  arising  from  the  virtual 
photon  processes  are  only  partially  cancelled  by  that  from  the  real 
photon  emission  process.  This  is  due  to  the  phase  space  restrictions 
imposed  on  real  photon  emission.  Therefore,  the  corrections  would  con¬ 
tain  a  logarithmic  term  which  diverges  at  the  high  energy  end  of  the 
pion  or  muon  spectrum.  This  term  is  expected  to  make  a  significant 
contribution  to  the  spectrum  of  pion  and  muon,  if  the  charged  lepton 
in  the  final  state  is  relativistic.  The  ratios  of  the  spectrum  with 
m^O  to  that  with  m^srO  for  both  muon  and  pion,  and  the  ratio  of  energy 
correlations  between  them  are  calculated.  A  rough  estimate  of  the 
effect  of  the  corrections  due  to  the  divergent terms  on  the  correspond¬ 
ing  ratios  are  given  so  as  to  illustrate  the  effect  of  the  correction 
to  the  phase  space  ratio  for  determining  the  mass  of  the  neutretto. 
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NOTATION 


The  following  notation  will  be  used  throughout  this  thesis. 

1.  The  matric  g^yhas  the  following  components. 

§>oo  ~  1  ^kk  =  — ^  for  kr  1,  2(  3* 

gMy  =  0  for 

2.  The  relativistic  notations  are  as  follows: 

(y0)2  =  -(rk)Z  =  1  k  =  l,  2,  3. 

r5  =  Y0rir2r3 

3.  Natural  units  have  been  used  i.e.fi=c  =  l. 

4.  The  notation  for  the  4-momentum,  3-momentum,  energy  and  mass  of  the 
particles  is  as  follows: 

4-momentum  3-momentum 

0 


(y5)  =  -i 


K 

P 

k 

q 

k« 


P 

k 

q 

k‘ 


energy 

mK 

£ 

w 

£ 

w' 


mass 
mK 
m, 


M 

"St 

\ 


K  meson 
muon 
pion 

neutretto 

real 
photon 

virtual  +  r  + 

photon  ° 

5.  U  is  the  spinor  of  muon. 

V  is  the  spinor  of  neutretto. 

6.  G'  is  the  coupling  constant, 

7.  G  is  the  4-momentum  imparted  to  pion  and  neutretto,  G=K-p=k+q. 


.. 


8.  /\  is  the  cut  off  mass  used  in  Feynman  regulator. 

9.  A.  is  the  fictitious  photon  mass,  which  is  introduced  to  avert  the 
infra-red  divergence. 

10.  6  is  the  polarization  vector  of  the  real  photon. 

11.  a  =  e  /4tt  is  the  fine  structure  constant  =  1/137. 

12.  0  is  the  angle  between  p  and  k*. 

13.  01  is  the  angle  between  q  and  k‘. 

14.  ^  is  the  angle  between  p  and  q. 

15 .  £  is  the  angle  between  k  and  k 1 . 

16.  7]  is  the  angle  between  p  and  k. 
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CHAPTER  I  INTRODUCTION 


An  upper  limit  of  about  3  Mev  has  so  far  been  placed  on  the 
neutretto  (muon  neutrino)  mass^.  It  is  worthwhile  to  study,  theoret¬ 
ically,  a  process  that  would  give  a  much  lower  bound  on  the  neutretto 
mass  and  would  at  the  same  time  be  experimentally  feasible. 


(2345) 

It  has  been  suggested  by  several  authors  that  the 

upper  limit  on  the  neutretto  mass  can  best  be  obtained  from  an  obser¬ 
vation  of  the  shape  of  the  spectrum  of  one  of  the  final  state  particles 
in  the  three-particle  decays  of  the  type 


y  B  +  C  +  V 


near  the  high  energy  end. 

In  general,  the  spectrum  of  the  particle  dPCE.m^)  of  one  of 
the  final  state  particles,  B  or  C  with  finite  neutretto  mass  rn^O, 
would  be  different  from  the  spectrum  dP(E,0)  with  m^=0,  for  there  is 
less  phase  space  available  than  if  m^O. 

Let  us  calculate  the  phase  space  integration  in  the  leptonic 

+  o  + 
decay  K—  — >  tt  +  p~~  +  V 

A  aV-p-k-q) 

2E  2^  26 


(1) 


■ 
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Where  K,  p,  k  and  q  are  the  four  momenta  of  K  meson,  muon,  pion  and  neu- 
tretto  respectively,  and  E,  d  and  6  are  the  energy  of  muon,  pion  and 
neutretto  respectively.  If  we  do  not  carry  out  the  p  integration  then 
we  would  get  the  muon  spectrum  dP(E,m^) 


dPCE.m^) 


f  3  3 

i  <Tk  <Tq 

j  2o)  2€ 


S^K-p-k-q) 


(2) 


To  work  out  the  integral,  we  set  G=K-p  (i.e.  the  four  momentum  imparted 
to  the  neutretto  and  pion  during  the  decay)  and  work  in  the  time-like- 
G-frame  defined  by  G=(GQ, 0,0,0)  (Appendix  A),  we  have 


dPtE.m^)  =  iLiG4 


„  2,  2  2v  ,2  2v2  2 

2G  (m^  +  m^)  +  (n^-m^) 


(3) 


In  the  case  where  m^=0,  if  one  evaluates  the  integral  in  the  frame  where 
K  meson  decays  at  rest,  i.e.  Ks(m^, 0,0,0) ,  then 


(k  -  p)2 


2  2 
mK  +  - 


2m^E 


(4) 


Since  we  are  interested  in  the  behaviour  of  muon  energy  spectrum  near 

2 

the  end-point,  G  should  be  close  to  its  minimum  value.  This  follows 

o 

from  eq.  (4)  where  G  is  minimum  when  E  is  maximum.  The  minimum  value 
2 

of  G  can  be  evaluated  from 


2  ,2  .2 

G  =  ( K  -  p)  =  (q  +  k) 


(5) 


,  . 
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The  right  hand  side  of  eq.  (5)  has  a  minimum  value  m  for  m.,=:0.  If  E 

tt  v 

2 

is  slightly  less  than  the  maximum  value,  say  E  =  E  -6E,  then  G  will 

lUcLJv 

2  2  2 

be  slightly  above  ,  say  G  =  G^  +  a;  therefore  dP(E,0)  close  to 
the  Emax  will  behave  like 


dP(E,0) 


tt  [(  2  x 2  2  x  2  41 2 

TT  ^(mtr  +  A)  -  2(mTT  +  A)mTT  +  “ttJ 
mTT 


TT 


2  i 

(a  ) 


(6) 


i.e.  proportional  to  the  first  power  of  A. 
If  we  had  not  set  my=0,  then 


TT 


dPCE.m^)  »  ? 

24 


",  2  ax2  2  Aw  2  2.  ,2  2x2 

(m^  +  A)  -  2(mTT  +  A)(mTT  +  my)  +  (m^-m^) 


JL  (A2  -  4m2m?.)2  (7) 

2m2  17  v 

TT 


In  the  above  expression  terms  of  higher  order  in  m^  have  been  neglected. 
From  eq.  (6)  and  eq.  (?)  we  can  get  the  phase  space  ratio 


dP(E,mv) 

dP(E,0) 


4m2  m2 
tt  V 


(8) 


This  can  also  be  expressed  in  terms  of  6E2  instead  of  A2.  From  eq.  (4) 


4 


JZ  ^  _2  r2 
dIta  •+■  m,.  -  (i 

e  =  _£ _ t _ 

2m^ 


(9) 


2  2 

Now,  when  G  =  G  +  A,  then 

min 


E  = 


mK  - 


2  2 

mn  +  nV 


2m 


K 


~  -  -  4E 
ZaK 


(10) 


therefore 


and 


6E  =  Z\./2m^ 

dP(E,m^) 
dP(E,  0) 


1  - 


2  2 
m  m.> 
tt  V 


m^(6E)2 


(U) 


(12) 


This  expression  is  quite  general  for  all  the  three-particle  decays. 

From  eq,  (12)  we  see  that  the  muon  spectrum  cuts  out  at  a  distance 
6E  =  m^^/m^  from  the  end  point  for  the  zero  mass  case.  The  linear  de¬ 
pendence  of  6E  on  m^  is  an  experimental  advantage.  On  the  other  hand, 
the  dPCEjm^)  differs  from  dP(E,0)  also  due  to  the  dependence  of  the 
matrix  elements  on  m^,  if  m  ^0.  But  if  it  is  assumed  that  the  leptons 

couple  locally  and  the  leptonic  current  is  of  the  type  UYu(l-iYc)V,  then 

r  j 

the  contribution  to  the  spectrum  arising  from  the  ra^-dependent  term  in 
the  matrix  element  is  always  quadratic  in  m^  while  the  contribution 
arising  from  the  phase  space  considerations  alone  is  linear  in  m^.  Thus 
for  small  m^  the  phase  space  contribution  dominates  the  contribution 
coming  from  the  matrix  element. 


■  ■  I  1  X..  ' 

•  c.-'  ‘  OO  5  \i 

- 

*  " 

*'  ;  - 


5 


If  one  plots  dP(E,iO/dP(E,0)  versus  energy  of  the  observed 
particle,  and  looks  for  the  deviations  of  the  spectrum  at  the  end-point 
from  the  horizontal  straight  line  appropriate  to  m^srO,  a  limit  on  m^ 
can  be  set* 


Six  different  spectra  in  the  decays  A0  — >  p  +  jj~  +  y, 

+  V  +  y,  K"  — >  tt°  +  }JT  +  V  were  suggested  in  both  the  papers  by 
Dennery  and  Primakoff  (1963)^  and  Ginsberg  (1965)^.  Greatest  devia¬ 
tions  may  be  anticipated  in  the  spectra  from  the  first  two  decays,  but 

because  of  the  relatively  short  lifetime  and  electrical  neutrality  of 
o 

the  A  and  because  of  the  difficulty  of  making  precise  measurement  of  y 

+  O  “I" 

energy,  the  spectra  from  rr  +  p”  +  2;  would  be  more  favorable  to 

study  experimentally* 


To  see  the  effect  on  the  spectrum  of  radiative  correction  on 
the  decay,  let  us  take  the  case  where  there  is  a  real  photon  being  emit¬ 
ted  with  four  momentum  k'  (innerbremsstrahlung).  If  E  is  close  to  E  , 

IUaX 

that  is  to  say  we  only  look  at  the  end-point  energy  in  the  spectrum  at 
which  the  effect  on  is  significant,  then  the  photon  must  not  be  al¬ 
lowed  to  carry  away  a  lot  of  energy,  and  must  therefore  be  a  soft 
photon*  The  phase  space  integral  then  becomes 


3  3  3  3 

dp  d^k  d  q  d  k1  c4/t,  .  ,  ,  x 

— £. - 2, -  6  ( K-k’ -p-k-q) 


■j 


2E  2a)  26  Zoo7 


(13) 
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Following  the  same  procedure  as  used  before,  with  the  replacement  of  G 
by  G'  =  G-k1 ,  we  have 


IT 


dP(E,m  ,)  =  — - 


2G‘ 


2 


d^k' 

2co' 


(G1^  -  2G,2(m2+m2)  +  (m2-m2)2 


1 

T 


(14) 


To  illustrate  the  importance  of  the  radiative  corrections.  It  is  suf¬ 
ficient  to  deal  with  the  case  my=0.  One  has  by  energy-momentum  conser¬ 
vation 

(G  -  k' )2  =  (q  +  k)2  (15) 


In  the  time-like  G-frame  G  =  (GQ ,0,0,0),  eq.  (15)  becomes 


Go  -  2kiGo  =  ^  +  k)2 


(16) 


(k^  is  maximum  when  (q  +  k)  is  a  minimum  or  G  is  a  minimum.)  If  we 

2  2  t 

allow  G  to  be  rn  +  A  as  done  earlier  the  maximum  of  k^  is  determined 

TT  O 


by  eq.  (16)  as 


%  +  A  -  2ko 


m 


G  =  m2 

O  TT 


or 


=  A /2m 

max  tt 


(17) 


A/  2m 


dP(E, 0) 


rr 


k' 


(G-k«r  - 


2(G  -  k'  )2m2  + 


i 

ip* 


K 


(18) 


by  expanding  the  terms  in  the  bracket,  and  neglecting  terms  of  order 

1 2  2  2  2 
k*  zv  in  favour  of  terms  of  order  kT  or  ,  we  get 
o  o 


7 


dP(E.O) 


A 

2m 


TT 


k» 

o 


A.  +  AG  k‘ 

o  o 


4aG  k« 
o  o 


dk» 

o 


A 


p2ni 

J  k'  (Term  of  order  A)  dk’ 


A" 


Thus  in  this  case  dP(E,0)  is  proportional  to  the  third  power  of  A. 

This  result  is  quite  different  from  the  case  where  the  real  photon  is 
absent.  Because  of  the  lack  of  phase  space  available  to  the  real 
photon,  the  emission  of  real  photon  is  inhibited  at  the  end  point. 
Therefore  the  infrared  divergent  terms  arising  from  the  radiative  cor¬ 
rections  of  the  virtual  photon  processes  will  be  only  partially  can¬ 
celled  by  those  from  the  real  photon  emission.  Thus  the  corrections 
would  reflect  this  lack  of  cancellation  of  the  divergent  terms  through 

an  infinity  at  the  end  point  of  the  spectrum.  It  has  been  shown  (Chap- 

2  mm 

ter  V)  that  such  a  term  does  indeed  appear  in  the  form  —  B(3)log(l - JLJijt 

I  4tt  m£6E 

where  B(3)  =  l+^log—-^  with  3=p/E  which  is  finite  for  the  range 

of  momenta  available  to  the  muon,  while  log(l-  acs)  diverges  logarith¬ 
ms  6E 

mically  at  the  end  point  of  the  muon  spectrum.  This  effect  was  first 

anticipated  by  Allcock  (1965)^  his  w0**k  on  radiative  corrections  to 

(4) 

the  decay  jj\  ■ — ^  e  +  v  +  Later,  Kamal  (1965)  worked  on  the  decay 
tt  — >  /J  +  V  +  y  and  showed  that  this  correction  term  will  not  have  an 
appreciable  effect  on  the  ratio  dP(E,m^)/dP(E, 0) .  This  happens  because 


■ 


£  f'-  -*<" 

■ 
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in  the  tt^  case,  the  charged  lepton  is  not  relativistic,  or  in  other 

words,  the  radiative  corrections  are  damped  out  by  the  small  value  of 

,  l-fc^ 

function  B(3m)  =  1  +  — =—  log  - — - — .  He  also  suggested  in  the  paper  that 
the  radiative  corrections  would  be  expected  to  be  large  in  the  case  of 
K  ^  decay,  for  the  mass  of  K  meson  is  large,  and  therefore  the  charged 
lepton  can  become  relativistic  (i*e.  3^  1),  Hence  B(3m)  would  be  of 

order  unity.  Consequently  the  logarithmically  divergent  term  will  not 
be  damped  out  and  would  be  important. 


In  Chapter  II  a  brief  account  is  given  of  the  derivation  of  a 
model  interaction  Hamiltonian,  The  model  Hamiltonian  is  justified 
through  the  use  of  the  intermediate  vector  boson  hypothesis.  With  this 
assumption  some  simplifications  are  made  for  the  convenience  of  later 
calculations ,  Our  results  are  expected  to  be  independent  of  the  form 
of  the  Hamiltonian. 

All  virtual  photon  corrections  are  made  in  Chapter  III  and 
real  photon  corrections  in  Chapter  IV, 

The  ratios  of  the  spectrum  with  m^O  to  that  with  m^=0  for 
both  the  muon  and  pion  are  calculated  in  Chapter  V,  The  energy  cor¬ 
relation  between  the  pion  and  muon  is  also  dealt  with  in  this  chapter. 
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CHAPTER  II  FORM  OF  THE  AMPLITUDE 

The  decay  mode  of  K-— +  tt°  can  be  described  by  a  Feyn¬ 
man  amplitude  whose  general  form  is  restricted  by  Lorentz  Invariance 
and  Time-reversal  Invariance.  We  shall  assign  a  local  leptonic  current 
of  the  (V-A)  type,  U^y^l  -  iy^)^,  the  l©pton  pair.  The  decay  pro¬ 
cess  then  looks  like  fig.  1. 


The  strong  interaction  involved  in  the  transition  to  the  intermediate 
state  with  the  emission  of  the  pion  is  not  known.  But  we  can  write 
down  phenomenologically  a  current- current  type  interaction.  The  decay 
matrix  is  then 


<K |  gJa | TT°>U^ya(l  -  ir5 ) V  =  gfc(K  +  k)a  +  d(K-k)Juya(l-ir5)V  (19) 


where  c  and  d  are  functions  of  G  =  K-k,  and  are  real  if  time  reversal 
invariance  holds  for  the  decay  interaction;  g  is  the  coupling  constant. 

We  can  make  some  simplifications  to  the  Feynman  amplitude  by 


1 


■ 
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the  assumption  that  an  intermediate  vector  boson  state  exists 


the  diagram  would  be 


2 


Then, 


Assuming  a  phenomenological  coupling  between  the  boson  B  and  the  (rrK) 
system,  the  contribution  of  Fig.  2  to  the  decay  matrix  would  be  of  the 
form 


(K  -  k) 


Drtt(i-ir5)v 

a  (K-k)2  .  m2 
B 


2  2 

If  nig  »  (K  -  k)  then  the  denominator  is  essentially  a  constant  and 
the  decay  matrix  can  be  simulated  by  a  phenomenological  weak  inter¬ 
action^ 

g,(Wk- W^(1-V  •  (20) 


where  G'  is  the  coupling  constant.  Note  that  equation  (20)  is  analo¬ 
gous  to  equation  (19)  if  the  constant  d  in  equation  (19)  is  zero. 


The  electromagnetic  interaction  is  introduced  through  the 
replacement  of  3 ^  by  +  ieA^  in  equation  (20).  Hence,  we  have 


,  ,  '  ’’  f 


■ 


. 


11 


G' 


^ttH1  +leA/u^K"^Kr)j  ^rr]  %i  ^ 1 "  ^j) v 


(21) 


In  equation  (21)  the  operator  operating  on  ^  is  unchanged,  because  the 
pion  is  electrically  neutral.  By  expanding  eq.  (21)  we  see  that  there 
is  an  additional  term,  which  is  absent  in  eq.  (20). 


ieG'*!ITA?J«iKUY/ja-ir5)V  (22) 

This  term  corresponds  to  the  diagram 


Fig.  3 

TT 

i) 

As  a  consequence  of  this  result,  the  diagrams  of  the  virtual 
photon  processes  with  the  photon  emitted  at  the  vertex,  like  fig.  5 
and  fig.  7,  should  be  evaluated  and  included  in  the  corrections. 


It  is  expected  that  the  results  we  have  on  the  ratio 
dP(E,ml/)/dP(E,0)  would  be  independent  of  the  model.  This  expectation 
can  be  made  plausible  on  the  grounds  that  had  we  chosen  the  more  general 
matrix  element  of  eq.  (19)  the  important  radiative  corrections  would 
appear  as  a  multiplicative  factor  to  the  entire  matrix  element,  leaving 
the  ratio  dP(E,m^)/dP(E,0)  unaltered.  The  important  radiative  correc- 


'  ^ 

. 

it  .r 


. 

/■ 

i  ; 
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tions  come  from  the  innerbremsstrahlung  graph  where  the  photon  is  emit¬ 
ted  from  the  charged  lepton  leg  (see  Chapter  IV). 


.  •  ■■  ■  .  ■ 
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CHAPTER  III  VIRTUAL  PHCTON  CORRECTIONS 


The  perturbation  theory  will  be  used  in  the  calculation  of 

the  radiative  corrections  of  the  K1 - ?  j-l"  +  tt  +  V  decay.  We  are 

interested  only  in  the  radiative  corrections  of  order  e  ,  because  con¬ 
tributions  from  the  higher  order  terms  are  small  enough  to  be  neg¬ 
lected  (Allcock  (1965) 

3.1  THE  SELF-ENERGY  OF  K  MESON 


First,  let  us  consider  the  self-energy  of  the  K  meson.  The 
Feynman  diagrams  corresponding  to  this  process  are  fig.  4  and  fig.  5. 
The  reason  for  inclusion  of  fig.  5  has  been  given  in  Chapter  II. 


The  matrix  element  corresponding  to  fig.  4  is: 


-iG'UY/J(l-iY5)(K"  +k)' 


Ml  1 


(-le)(K"+K')P 


1  1 
(2*)^  K'Z  -  m2 


(2n)4  (K"2-m|) 


•  ■  ■ 


■  :  • .  ■ ....  !oo;a  '  J  ‘  -  j  ...  ft.  a 
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~  (-ie)(K«  +K)pV(2rr)46(K»'-p-k-q)(2TT)46(K'+t-K)(2TT)4 
(2tt)^  \T  -  r  r 

6(K-t-K' )d4Knd4K‘d4t.  (23) 

In  the  above  expression  we  have  omitted  the  normalisation  factors  such 

i  /5  ii 

as  <5/0  /  F”  for  the  fermions  and  “Tm  s  for  the  bosons, 

(2tt)3/Z  y  i  (2tt)3/Z  yiij 

As  we  are  going  to  take  the  ratio  of  the  spectrum, these  factors  even¬ 
tually  will  be  cancelled.  After  the  trivial  integrations  are  done  in 
eq.  (23)  •  we  have 


4 

V6(K-p-k-q)d  t  (24) 


The  difficulty  with  this  expression  is  that  we  have  K^-  in  the  de¬ 
nominator  and  this  is  zero  for  a  free  particle.  To  overcome  this  dif- 

(7) 

ficulty  we  introduce  a  damping  factor  g(t)  ,  replace  e  by 


r°  -iflt 

eg(t)  s  e]G(ro)e  dPQ,  set 


r  =  (r°, 0,0,0) 


and  generalize  by  eg(t)— *e  G(r)e“^^dr 


OO 

with  the  normalization  g  (0)  =  jG(ro)dro  =  1 

_oo 


.. 
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To  do  the  complete  calculation  let  us  write  down  the  T  matrix 
explicitly. 


(2tt) 


_  -ipxo  H  ,  -iK"(x?-x,J  -ikxo 

(-iG')Ue  ^  Yu(l-iYt:)(K',+kr  — e  3  e 

r  5  K»  -mj 


n  -iK‘(x, -x^)  -it(x-.  -  x~) 

(-ieg^))  2  e  1  ^  e  ^  (-ieg(t2)) 


K'  -  ra‘ 


K 


t  -  V 


iKxl  P  -i(ix~  444444 

e  (K"  +  K'  ^  (K1  +  K>  Ve  J  d  K"d  K'd  td  x^  Xgd  x^ 


(2tt) 
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-iK"(x2-x3)  -it^-xg)  -iK'Cxi-xg) 
J  e  e 


i2  .  ^ 


T 


K 


-irlX  -iT £X 


GC^e  G(r2)(K'«+K')p(K'  +K)rV 


p  -i( p+k+q-K)x^  ^ 


e 


d  K" 


4  4  4  4  4 

d  K’d  td  x-jd  Xgd  x^dr^dr^ 


Carrying  out  all  the  trivial  integrations,  we  have 


-e2G‘ 


%(i-ir5)(K+k-  Ti-  r2)HG(r1)G(r2)dr1dr2  ni  pox2  _2 


(K-ri-P2r -m‘ 


K 


- ^-5 - r  d\(2K-  t  -  2P,  -  r2)0(2K-  t  -  P  fsCK-P^-p-k-q) 

2-X2  (K-t-/^)2-  P 


(25) 


,  •  -nr  '  '  »  so  '  ■  ' 


«sh  -  £  '  — 
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Note  that  if  p^=p2=0.  We  get  back  to  the  result  of  eq.  (24).  We  may 
write  eq.  (25)  as 


-iG«  |  uy^(i-iY5)(K+k-r1-r2y"dr1dr2G(r1)G(r2) 


2(r , K) 64( K-ri-r2-p-k-q)V 


(26) 


Where  2  is  the  t  dependent  term 


£  (T.K)  =  -  -^-T  d4t  -1  ?  - - — s — t  (2K-t-2r‘r  )  (2K-t-r,)P 

(2tr)  j  t2-  kz  (k-t-^r-m2  1  2  P  1 

We  expand  2  in  powers  of  P  and  retain  only  terms  of  order  f. 


2  (P.K)  =  - 


ie< 


(2tt) 


d\  1 


t2-  \2  t2  -  2Kt 


- -  (2K-t)2+ 

(2tt) 


4  1 

d  x 


t*-l2 


--i - (2K-t)(3r!+r2)  -  d\ 

t2  -  2Kt  ^  (2tt r  1  t2  -  X2 


2ie2  Pa  1  (2K-t)2(K-t)  r/ 


(t2  -2Kt)2 


=  ^  +  x3 


(27) 


The  first  integral  1-^  is  the  self-energy  integral  independent  of  T  and 
is  a  function  of  K  alone.  This  function  must  be  a  function  of  K2  so 
that  it  be  Lorentz  Invariant.  But  K2  =  m^2,  hence  the  1^  is  a  constant 
term  and  would  be  an  infinite  constant.  It  is  a  quadratically  divergent 
term,  which  can  always  be  cancelled  by  a  term  -6m  . 


a 


. 


■ 
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Therefore  the  integrals  left  for  evaluation  are  and  I y  By  using 

/  q} 

the  Feynman  cutoff  method'  ,  we  get 

Iz  =  ■  i?  log  i + $  v3^ +rz)p  (28) 

K 

j.  ,2  4  p 

I,  =  -2-J  (log  ^  -  2  log  —  +  6)2K  ry  (29) 

16tt  m^  X. 

Substitute  these  back  into  eq.  (26)  and  eq.  (27) 


-iG' 


J 


UY„(l-iY5)(K+k-rrr2)PV  dridr2G(r1)G(r2)  - 

Oc-fj-rg)  - 


m 


K 


I^[('  ?  106  5  '  I)Kp(3ri+r2)P  +2(log  j  -  21°®  ?  +  6)Vi P] 


K 


K 


(2tt)  o(  K-r1-r2-k-P-q) 


(30) 


We  now  symmetrize  the  factor  P  in  eq.  (17)  by  writing  — t>-§T^  +  p^. 

Then  eq.  (30)  becomes 
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•iG-  DYf/i-iY5)(K+k-r1-r2)''v  dr1dr2G(r1)G(r2)  -(r  ~ , 

s/’  X  fc 


W !  1‘ I  ^  I  ‘  I,2Vr.-"i>F*  - 

K  K 


A 


m 


2  log  -*-  +  6) 

r 


K  (Pi+r2)^>  (2TT)/l6(k-ri-r2-P-k-q) 


(3D 


In  the  limit  of  r^,  — >0,  we  get 


Tx  =  -  iG'UYH(l-iY5)(K  +  k) 


M 


-  <loe  4  *  t> 

4n  A.2  4 


(2n)  6(K-p-k-q)  (32) 


e  1 

Where  - —  =  a  is  the  fine  structure  constant  «  — — 

4tt  137 

A  is  the  cut  off  mass  used  in  Feynman  regulator 
A  is  the  fictitious  photon  mass,  which  is  introduced  to 
avert  the  intra-red  divergence.  As  A— >0  the  tern  log  A 
is  infinite. 

Now,  we  consider  the  correction  due  to  the  process  of  fig.  5. 
The  matrix  element  corresponding  to  this  process  is 


u 


(a^^CK'+t-p-k-qKa^^^SCK-t-k*  )d4K‘d4t  . 


' 


■ 


-  - 
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After  carrying  out  the  K’  integration  with  the  help  of  the  delta-function, 
we  get 


•e2G' 


Uya(l-iy5)  - — - j— - -  (2K-t)aV6(K-p-k-q)d>t 

( K-t )  -  m^  t  -  A 


(33) 


The  above  integral  will  not  give  rise  to  an  infra-red  divergent  terra. 
By  using  the  Feynman  cutoff  method,  we  get  the  result: 


2 

G'U(|  log  f—  +  2)|t(l.iY5)V(2TT)46(K.p-k-q)  (34) 

mK 

where  A  and  a  are  the  same  quantities  as  defined  earlier. 


3.2  THE  SELF-ErJBRGY  OF  THE  MUON 


Similar  to  the  case  of  K  meson  self-energy,  we  have  two  types 
of  muon  self-energy  graph,  fig.  6  and  fig.  7. 


Fig.  6 


■ 


■  -  - 


'  ,  .1 
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The  matrix  element  corresponding  to  the  fig.  6  is: 


t2  = 


(-1G'  )U(-ier  ) 


(7\  l 


1  ,  ,  x  i  1  -i  1 

(-leyJ 


(2tt)4  P'  -m/>  a  (2rr)4  ^-"V  (2 n)4  t2  -  X2 


ruU-irc)(K44t)4V(2TT)46(K-p"-k-q)(2TT)46(p"-t-p')(2TT)4 


/  v  4  4  4 

6(p'+t-p)d  p"d  p'd  t  . 


(35) 


After  carrying  out  the  trivial  integrations. 


t2  = 


'  p  _  a 

-e  G'Uy 


^  m 


A1  1 


(p-t)2-  m2  a  p2-  m2  t2  -  A2  ^  ^ 


y  (l-iy.XK+k^V 


6(K-p-k-q)d  t  . 


(36) 


2  2 

For  a  free  muon  p  =  m^.  Therefore  the  denominator  is  identically  equal 
to  zero.  We  may  handle  this  in  the  same  way  as  before,  i.e.  by  re- 


oo  ip  t 

placing  e  by  eg(t)=e  J  G(p)e  dp.  Then  the  matrix  element  becomes 


«0O 


.f  V.v! 


21 


(2tt) 


12 


_  -iPxl  a  1  -ip'^-x^  x  -ip’Kx^-Xg) 


e  G'U  e  y 


t- 


m, 


-  Mi 


-irixl  “ir2x?  l  -itCxg-Xn)  w 

G(rx)e  G(r£)e  —  ?  e  1  Y/J (1  -  iy5)( K  +  kf  V 


t2-A2 


i(K-k-q)xo  IlIiLlLiLlLlULl 
a  x^d  X£d  x^d  p"d  p‘d  td  P^d  l£ 


=  -iG» 


uy<P.D  .....  .  I  -  Yll(l-iY5)(K+kf  V(2n)46(K-p-k-q) 


G(r1)G(r2)dr:Ldr2 


(3?) 


Where  (p,D  is  the  t  dependent  integral, 


(2tt) 


E(p^)=7=V  ^y0— L — r. 


J 


jL^+f-j+m  0  t2-A.2 


This  is  a  logarithmically  divergent  term.  As  before,  we  expand  ^~|in 
powers  of  T  and  retain  terms  of  order  P. 


V(p.r)  =  - 


ie  ,4.  a  1  1  i©  ,/k  a  1 

d  ty  -  y  _  — -  + - -  d  *tyu  - 

°  t2-\2  (2tt)^ 


(2tt) 


1  1  n 

/  y  Y<?  .2  ,2  ^ 
p-t-m^  t  -A. 
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=  Ii*  +  V 

t 

The  term  1^  corresponds  to  the  diagram  £_  p 

p-t 


■7- 


(38) 

which 


would  be  cancelled  by  the  term  -  6m^.  The  integral  1^  can  be  evaluated 
by  standard  methods  as  was  done  before. 


s  -  £  t®-  i  •  ^  ■ ■  "•*$  -  if  - 1> 


(39) 


In  which  A  ,  A.»  a  are  the  same  quantities  as  defined  earlier.  Substitut¬ 
ing  the  expression  back  into  eq.  (37)»  we  obtain 


t2  = 


.iG1 


4 +  ¥  - 

L  y 


m 


-  log 


M 


>l 


- I - V  (l-ir5)U+k)HV(2n)  6(K-p-k-q)G(r1)G(r2)dr1dr2 


(40) 


We  now  symmetrize  the  factor  P  in  eq.  (40)  by  writing  P  =  2'(Pi+P2^  an<^ 
recalling  that  ^-m^  operating  on  the  free  particle  Dirac  spinor  gives  zero, 
we  may  add  the  operator  j(  -  m^,)  to  the  first  term  in  the  bracket  in 
eq.  (40)o 

.2 


-iG 


— /  cl  v 
U(-  2^ 


Z  —g  +  J-  2  (l^+/] ]_+/" 

m 


( ^+Pi+/2"mp|) 


rH(l-r5)(K+lc)HV(2iT)'''6(K-p-k-q)G(r1)G(r2)dr1dP2 


(41) 
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As  rif  r 2  ■ — >  o,  the  above  expression  becomes 


-iG'U 


-  (flog  A  +  5) 


4tt 


m2  4 

fA 


Yf)  (l-ir5)(K+k)MV(2rr)45(  K-p-k-q) 


('12) 


We  rewrite  the  second  term  as  follows: 


77/  a 


-iG«  U(£-) 
2rr 


2 


m 


,^(iog  3  - » 


^+/1l+/12+mH  . 


- 7P.(rl+r2)Yu(l-iYc- 

2p.(  1+2)  1  2  ^  3 


(K+kf  V(2rr)46(  K-p-k-q) 


In  the  limit  P^,r2  — >0,  this  reduces  to 


-iG'U 


a  / ,  “J  ,  x 
^(log  -£  -  1) 


Y^j  (l-lr5)(K+k)HV(2rr)  ^(K-p-k-q) 


(^3) 


Adding  eq.  (42)  and  eq.  (43),  we  have 


T  - 
2  ~ 


-iG'U 


a 

4tt 


m 


M  _  lino-  __  _  2. 


(log  -4  -  flog 


2 

“jji 


f) 

4 


y  (  1-1y5)  (  K+k)HV(2TT)46(  k-p-k-q) 

(44) 


The  matrix  element  for  the  graph  7  is  given  as  follows: 


U(-ieya) — (-ieG'y^)  “1 


(2tt)4  i- 


(2TT)if  t2-X2 


(i-ir5)v(2Tr) 


6(  k-p '  -q-k-t )  (2tt)  *5( p-p ' -t  )d\d^p ' 


(45) 
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f  „  0  W+m 

=  ~eG!  \Uy 


(p^t)2-!?!2  °  t2~\‘ 


1  (l-iy^)V6(K-p-k-q)d4t 


The  above  integral  does  not  contribute  any  infrared  divergent  term. 
After  integrating  with  the  help  of  Feynman  method,  we  get 


~  G3U(31og 


A  ^  O  If, 

+  — )m^(l~iy^)V(2Tr)  6(K-p-k~q) 
m^ 


(46) 


3.3  THE  VERTEX  CORRECTION 


The  Feynman  diagram  corresponding  to  the  vertex  correction 
is  given  in  fig.  8. 


The  expression  for  the  matrix  element  for  this  process  is  as  follows s 


(-iG!)U(-ieya) 


J 


(2tt)4 


-i-  y  (l-iY5)(K>+k)M 


i  1 

(2n)lp  K*2-]*2 

A. 


(”ie) 
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U+K')  — -4  V(2TT)46(K‘-p»-k-q)(2TT)46(K»+t-K)(2rr)4 

a  (2rr)4  t2 


6(p,+t-p)d^psd^K,d\ 


=  -  e2G»Uyu  — i -  y  (1-iYcXK-t+k/  - i-w — x  (2K-t)  ~  (K-p-k-q)d\ 

J  i-t-^  M  3  (K-t  )2-m2  Qt2 


=  (~i)(2rr)Vui(K,p)V6(K-p-k-q)  (4?) 

Where  K,p)  is  the  p-dependent  integral  whose  explicit  form  is 

-  (2K-t)a  i  d\ 
mK  1 

(48) 


J(K,p)  = 


-1 


(2tt) 


Y 


i-U 


1  y^(k-t+k)M  (l-iy^) 


m, 


(K-tr- 


The  numerator  in  eqa  (48)  can  be  arranged  as 


(49) 

The  numerator  being  of  order  the  integral  _^(  K,  p)  is  logarithmically 
divergent.  Introducing  a  regulator  mass  A  and  carrying  out  the  standard 
Feynman  Parametrization  the  denominator  of  ^(K,p)  can  be  written  as 


r 


t2  t2~2pt  t2-2Kt 


d\  = 


A 


-6 


1 

r 


dL  dy 


KJ 

o 


u(l-u)dudt 


u 

o 


( t-kyu)2-kyU2-L(l-u) 


(5°) 
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where  k  =  p  +  G(l-y).  Substituting  both  eq.  (50)  and  eq.  (49)  into 

y 

eq.  (48),  we  have  the  expression, 


yu,P)  = 


-i 


(2tt)‘ 


J 

42 

1 

i 

r 

"\ 

f-6 

dL 

dy 

u(l-u)du 

k.  V. 

2 

\T  a 

vj 

n 

dt 


(  4k.  p-2  M-U-Un+iz  )(t+t) 


(t-k  u)2-k2u2-LCL-u)]‘ 


A 


+6 


1  1 

r  p 


dL 


d  y 


u(l-u)du 


A.  o  o 


dt 


4K.ptf-2$2-t!#-l!2iTiM+)(3 


vy 


(t-k  u)2-k^u2-L(l-u) 


.(i-ir5) 


On  Shifting  the  origin,  by  replacing  t  by  t+k  u,  and  dropping  all  terras 
of  odd  order  in  t  in  the  numerator,  we  obtain 


All 


^(K,p)  = 


(2tt  y 


-6 


dL 


dy 


u(l-u)du 


A  o  o 


dt 


( 4K.  p-J^u-ny,  j(yU-2  J^u+k^u2  ) 
|t2-k^u2-L(l-u)]  4 


All 


-6 


p 

r\ 

t2 

dL 

dy 

u(l-u)du 

dt  " 

X* 

o  c _ 

0 

o 

t2-k^u2-L(l-u) 

0W) 


A~  1  1 

f\  ^ 


+6 


dL 


dy 


r 


4K.  p^yu-2^u2-^y^yu^-ra/u^yu^+^u- 


u(l-u)du 


,2  J  ° 

A  o  o 


2  „  ,/2  2  w3  3 


dt 


t2-k2u2-L(l-u)l 

•X 


■ 
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All 

P  r>  r\ 


+6 


dL 


dy 


u(l-u)du 


dt  t 


-2  Kt2+j^t2 -m^t2 


A  o  o 


t2-k|u2-L(l-u)J 


All 


-l 


(2tt) 


-  |  dL  I  dy 


A  o  o 


u(l-u)du  , 

[(k|u2+L(l-u)] 


1  1 

2  s  (ky.p+^y)u 
-f2iTT  dL  dy  u(l-u)du 


A  o  o 


[k^u2+L(l-u)]: 


All 


-ITT 


dL 


dy 


u(l-u)du 


.2  2 

kyu 


~  O  ° 

A^  o  o 


Jk^u2+L(  l-u)j 


+2irr 


A 

1 

1 

1 

dL 

dy 

u(l-u)du  - 

p2  ^ 

V. 

k^ir+L(l-u) 

“  y 

(K+k) 


+1TT 


A2  1 

1 

^  p 

4K.pl/ u 

dL  |  dy 

u(l-u)du 

y 

r  p  o  v~i 

J 

J2 

[kV+L(l-u)J 

r  l  l 


-irr 


dL 


dy 


u(l-u)du 


A  o  o 


(  2  )  U2 

[k^u2+L(l-u)]  2 


(i-iy5) 


. 


4 
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All 

r\ 


+iTT 


dL 

kJ 


dy 


VJ 


A3 

u(l-u)du  y 


k  o  o 


[k?u2+L(l-u)J 


All 


-f2in 


dL 


dy 


u(l-u)du 


2j4im. 


"■’2  ^ 

AT  o  o 


[k^u2+L(l-u)] 


d-ir5) 


(51) 


In  eq.  (51)  only  the  first  integral  diverges  when  k= 0.  Therefore  for 
the  rest  of  the  integrals  we  may  take  0. 


^(K,p)  = 


ie2 

16tt2 


-2K.p 


2  1 
\c  r 

4  log  -f  dy  +  2 

y  * 


dy(ky.p+^y)  —  -  j  + 


dy  log  ~ 
k2 

y 


i  (fetf)  +  4K.P 


dy-| 

ky 


l 

r 


i. 

■2 


: 

dy-|  + 

4 


+  2(2j(+2m^) 


dy(  log^)-i 


J 

o 


•  d-ir5) 


y  J 


(52) 


In  the  above  expression,  there  are  five  integrals  that  need  to  be 
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evaluated. 


1 


kyady 


1 


j 


logkydy 


P  kyokyrdy 


dy 


To  evaluate  these  integrals  we  make  the  substitution  y=l-z,  then 


G2(z2 


+  2 


p.  G 


2 

+  ~2> 
Qt 


If  we  take  =  a  and  0?)^  -  (m^/G^)  =  b,  then 
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2  2 

k  =  G( z+a+b) (z+a-b) 

When  one  evaluates  all  those  integrals  in  terms  of  these  variables, 
one  gets 


J 

o 


1 

_ dz _ 

G2(  z+a+b  )  ( z+a-b ) 
o 


log 


(z+a+b) (z+a-b) 
LXZ 


G2 

=  F  log  ~  +  F 
°  XT 


t 

o 


(53) 


Where 


log( 


1+a-b 

1+a+b 


■)(££) 

a-b 


(5^) 


Fo  = 


2^ 


||[log(l+a-b)]2-  j  [log (a-b  j2-  j  [log (1+a+b)] 


+  |-[log(a+b)]2-  log  (1+a+b)  log  (1+a-b)  + 


+  log(a+b)log(a-b)  +  L(-  3-±a~Jp) 

2b 


+ 


4- 


log2b.log(iiS=^ 

a+b 


(55) 


The  L  functions  in  the  eq.  (55)  are  the  Spence  functions,  which  are 


L(t) 


t 


r\ 


vJ 


log(l-x) 


o 


defined  as 


' 


■ 


. 
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Near  the  end  point  of  the  muon  (or  pion)  spectrum,  the  term  designated 
as  F*  will  not  give  rise  to  a  divergence.  Ely  making  the  same  substitut- 

Q  i 

ing  as  that  for  JQ,  we  can  evaluate  the  rest  of  J's. 


Ji  = 


p  +G  z 

a 


- - — - -  dz  =  p/0  +  G/i 

G  (z+a+b)(z+a-bj 


Where  F-,  =  — ~ 

2bG2 


-(a-b)log(2±a=k)  +  (a+b)log(i±2±8) 
a-  b  a+b  J 


(56) 


Jo  = 


log 


G2(z+a+b)(z+a-b)j dz  =  logG^  + 


Where  F„  = 
2 


(a-b)log  —?Fb  +  (a+b)log  +  log(l+a+b)(l+a-b) 

a-b  a+b 


(57) 


J3  = 


PaP  +(GaP  +PaG  )z+GQGrz' 
G2  ( z+a+b  )  ( z+a-b  ) 


dz  55  PaPrFo+(GaPr+PaGr^Fl+GaGrF3 


Where  Fo  =  — -k- 

J  2bG2  L 


(a+b)2log(--+aj*b)  -  (a-b)log(-1-t.arb)  +  3b 
a+b  a-b 


(58) 


J,.  = 


kyX  dy  =  pa+|Gc 


(59) 


Substituting  these  results  back  into  eq.  (52).  we  have 
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4(K,p)  =  JL-  |[-2K.p(F0log  2^  +  F^)  +  2^^^)  +  20^+1$^) 


A 


-  i  +  log—  -  f2J 


(tf+lO  +  4K.p(|(Fo+^F1)  +  Jt  -  mH  + 


2 

+  iii3F0  +  i&i^F1  +  i^^F3  +  iF4  +  (4)t+inp)log  ^  - 


-  (^Um/J)F2 


-  d-ir5) 


(60) 


If  we  are  working  in  the  frame  where  K  meson  is  at  rest,  then 
K  =  (m^, 0,0,0),  and  U  =  m^U.  So  we  can  write  eq.  (60)  as  follows: 


4(k,p)  =  S-, 
16tt 


2  ( r  q2  2  2 

-2raKE(Folog  —  +  F'o)  +  2(mpF0+raKEF1-m‘cF1)  + 
\ 

+  4mkE(F0-F1)  +  2m^+  log  L  _  1  .  f2]  (M) 


imn  -  2V°K  +  H‘\E)  +  4“/<”1KE(Fo-F1)  +  ^Foml  *  m?Fl 


2  A  -a 

-  zFfijA mK  -  m^log  +  4m/>|F2  +  fm^ 


-2m^  +  4mJSF^  + 


+  mKEFj  -  m2  +  41og  2.  _  4F2]  tf  -H  1-  (1-ly,) 

G 


(61) 


Finally  the  matrix  element  for  the  vertex  case  is 


T3=  (-i)G'u  £ 


(2mKEFolog\  +  f|)(K+l0  +  f j  +  f3^  + 
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(l-iY5)V(2TT)46(K-p-k-q) 


(62) 


Where  f|  =  -2mKE(F0logG2+F^)  +  2(m*F0-mKEF1-m^F1)  +  ^E^-F-^  + 


+  2mK*  log^-F2-l_ 


(63) 


f2  =  4lymKE  *  |  >  *  2V*K  +  HmKE(Fo*Fl)  +  ^Fo“2  ‘  m/«Fl  * 


2  3 

2DVtmKF3  +  2^3 


A 


m^log  “2  +  4lVF 2  * 

G 


(64) 


f^  =  -2m^  +  4mKEF1  +  m^EF^  -  +‘  4log  ~  -  4F2 

Qt 


(65) 


In  eq.  (62)  the  log  A.  term  is  the  infrared  term  that  we  had  before.  The 

i  2 

function  FQf  FQ,  F^*  F2  and  F^  are  functions  of  G  throughout  the  phys- 

2  2  2 

ical  region  (m^-m^)  >,  G  >.  m^.  At  the  end-point,  G  tends  to  m^,  and 

these  functions  vary  slowly  as  the  end-point  of  the  muon  (or  pion) 
spectrum  is  approached.  It  is  justifiable  to  regard  the  functions  of 
Fl*  F2  an<^  F3  as  cons^an^  when  compared  to  the  divergent  terms. 


Adding  all  the  virtual  photon  corrections  to  the  uncorrected 
decay  amplitude  T  ,  one  has 


T  =  T  +  T-,  +  T-I  +  To  +  Tl  +  To 

vp  o  1  1  2  2  3 


(66) 


where  TQ  =  (-i)G' Uy^  ( l-iy^)  (  K+k)M 'i/(2rT)4o(  K-p-k-q) 


(6?) 


' 
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which  corresponds  to  the  diagram 


K 


Fig.  9 


Putting  the  eq.  (32) »  eq.  (34),  eq.  (44),  eq.  (46)  and  eq.  (62)  into 
eq.  (66),  we  have  the  expression  for  the  amplitude  corrected  by  the 
virtual  photon  processes. 


Tvp  =  (-DG-UY^l-ir^jl  +JL  j[-(2^mKEF0)logX2  - 
+  f £  +  f^KH  +  K^yQqa  ))  V ( 2rr) ^6 (  k-p-k-q ) 


(K+k) 


=  T  1  + 


a 

4tt 


-2-2mKEF0)logr-f1 


jgj  (_i)G'U(f2  +  f 3K  +  K)Dy  ( 1-iYj ) V ( 2tt) "*6 ( K-p-k-q )J  (68) 
where  we  have  written 

f^  =  f]_  +  2’logA^  +  ^log^  -  3  (69) 

.2 


f2  =  f2 


A  3  1  m. 


31og  ~  | 

m^ 


(70) 


■  . 
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-  ,i  3n  "  9 

f3  =  f3  ”  S’  log  “2  ”  4 

m, 


(71) 


In  the  expression  (68)  the  infrared  term  (2-2m^EF0)logX  will  be  can¬ 
celled  by  a  similar  infrared  divergent  term  arising  from  the  real  photon 
emission  process.  The  calculation  of  the  real  photon  emission  correc¬ 
tion  is  done  in  the  following  Chapter.  The  ultra-violet  terms  in  eq. 

A2  A2 

(68)  can  be  written  as  log  — —  or  log  —  ,  then  these  terms  will  not 

m| 

give  a  large  contribution  to  the  spectrum  for  any  plausible  value  of  A. 
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CHAPTER  IV  REAL  PHOTON  CORRECTION 
(INNER  BREMSSTRAHLUNG) 


There  are  two  types  of  innerbremsstrahlung  graph,  fig,  (10) 
and  fig,  (11),  which  will  contribute  to  the  correction.* 


The  matrix  element  corresponding  to  fig.  (10)  is 


*  The  innerbremsstrahlung  graph  of  fig.  (12)  will  not  five  any  contri¬ 
bution  to  the  correction,  as  we  work  in  the  frame  where  K  meson  is  at 
rest.  This  can  be  seen  as  follows:  The  vertex  factor  for  the  photon 
emitted  from  the  K  meson  line  will  be  (2K-k‘).£>  where  £  is  the  polari¬ 
zation  vector  of  the  real  photon  and  is  orthogonal  to  k'  i.e.  k'.£=0. 
Ne  can  choose  a  gauge  where  £  has  only  the  space  components.  Therefore 
in  the  time-like  K  frame  K.£=0.  Hence  the  vertex  factor  vanishes. 


Fig.  12 
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TjB  =  (-iG)U(-ierp) - ■£ 


(2tt)3/2  (rf-|f)-: 


1  YH(l-iY5)(K+k)f'v(2n)lt6(K-k'-p-k-q) 


m, 


(72) 


The  matrix  element  corresponding  to  fig.  (11)  is 


T'je  =  (-iG')U(eYH)  ■  P,;,  (l-iY  JV(2tT)46(K-p-k-q-k' ) 


(2tt) 


572 


(73) 


In  both  eq.  (72)  and  eq.  (73) »  k’  is  the  4-momentum  of  the  photon  and 
if  the  polarization  4-vector  of  the  photon,  k‘  and  £  are  orthogonal 
to  each  other.  If  one  is  working  in  the  frame  where  K  meson  is  at  rest, 
then  K.£=0. 

Adding  the  two  matrix  elements  together,  we  have 


TIB  =  TiB  +  T,iB=(-1G'^7^37i 


(2tt  )• 


£(^-tf‘-Hry,) 

(p-k)2-m2 


(K+K) 


(i-ir5)v 


(2 tt)  o(K-k*-p-k-q) 


(74) 


By  making  use  of  the  relations  u  i  =  mMu 

and  t  i>  =  2p.£  -  i>t 
we  can  simplify  eq.  (74)  as  follows: 


T  =  (-1G'  )U  - 2__ 

IB  (2tt)3/2 


p.k’ 


+  t  (l-iY5)V(2n)t6(K-p-k-q-k')  (75) 


The  decay  probability  due  to  the  photon  emission  alone  is  given  by 
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IB 


2  d3p  d3k  d3q  d3k» 
2E  2w  Z€  Zoo 1 


r 


J 


IB 


2  £2  A  £a  Al  + 


2E  2 


2  • 


mil  2  d3p  d3k  d3q  d3k* 
IB  I  2E  2w  2€  260' 


Vj 


2Re 


X 


IB 


d3p  d3k  d3q  d3k' 
2E  Zto  ze  2(0* 


(76) 


where  X^  is  the  cross-product  terra. 


2  n ,  2 

y  e  Or* 


~1  r 


TR  =  r  D(-£ih)(jtrt()(l-iY  )V  VU-iyO/U 

s  A.  (2tt)3  L  P.k'  3  JL  3  - 


(2rr)  6(K-kf-p-k-q) 

(77) 


Since  we  are  not  concerned  about  the  spin  of  the  neutretto  and  the 
polarization  of  the  photon,’ we  have  to  sum  Xjg  over  the  final  spin  states 
and  over  the  polarization  vectors  of  the  photon.  And  the  same  applies 


to 


T  * 

iIB 


and 


T" 


IB 


Let  us  calculate  the  first  term  in  the  eq.  (76).  We  can  write 


it  as 


T' 

2  d^2  d3k  d^a  d3k«  _ 

T 

2  d3p  d3k  d3q 

e2  Ip. £ 

2  d3k« 

IB 

2E  Zu)  ze  2co» 

V. 

0 

2E  Z6  ZcJ 

(2tt)3  j  p.k' 

2  2 

(78) 


Where  TQ  is  the  uneorrected  decay  amplitude. 


2  =  (2tt)L'2 


U(j(+lO(l-iY5)v|2a(K-p-k-q) 


(79) 
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(9) 

Summing  over  the  polarization  vectors  of  photon  ,  we  have 


I 


P.  t 


D2-(p.k,2/^'2)j  b2  (i-u2cos2e) 


X 


p.k* 


(Ec,)«-p.k«)2 


2  ? 
(l-UScosG)^ 


(80) 


In  where  U  =  —  corresponds  to  the  velocity  of  the  photon. 

=  p|/i£  ,  which  is  the  velocity  of  the  muon. 

9  is  the  angle  between  the  momentum  of  the  muon  and  that  of 
the  photon, 

?  3 

d-V  .  .  .... 


To  integrate 


u 


2 

P*£ 

P.k« 

2 

as  the  z-axis  direction  and  do  the  integration  in  the  spherical  coordin¬ 
ates  (Appendix  B)^^.  We  have 


P«£ 

2  d\» 

p.k’ 

2  2W» 

+  cO)) 


(81) 


with 


+i 


(log 


-1 


max) 

X 


(1-jdx)2 


(82) 


and  CO)  =  -4BO)(log2-2)  +  |  L(@)  -  L(-@) 


1 

0 


L(i±£)  -  L(i=2) 


-  *  |logfe|(log  +  2) 


(83) 
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In  eq.  (83)  the  L  functions  are  Spence  functions  as  defined  earlier  and 
BO)  is  defined  as 

B(0)  =  l  +  ilogM  (84) 

BO)  can  be  shown  to  be  identical  to  the  term  (l-m^EFo)  in  eq.  (68), 

If  we  write 


a  = 


p.G 


(hir/e)-  1  +  B2 


G2  (mK/E)2-2(m*/E)+l-.|32 


K' 


,  2  2 
b  =  /(EE)  -  2£d  = 


tBmK/E 


2  G^  (m^/E  )2-2(mK/E)+l-ji2 


then 


■mKEF0  =  -mKE 


2bG* 


=  -i-  log 

2B  1+B 


2  2 

Near  the  end-point  of  the  muon  spectrum,  G  =  m£j.  +  2mrr6E,  the  photon  is 
not  allowed  to  carry  away  too  large  an  amount  of  energy.  The  maximum 
value  of  k’  can  be  calculated  from  the  relation 

(G  -  k')2  =  (q  +  k)2  (85) 

The  maximum  photon  energy  occurs  when  (q+k)^  is  a  minimum,  i.e. 

(m  +  mw)  .  In  the  rest  frame  of  the  K  meson  G=  (m„-E,  p),  we  get 
TT  S'  A- 


■ 
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k* 

max 


mK6E  - 
-  E  -  ^cos0 


(86) 


Therefore  we  can  write  I  as 

o 


)%6E  -mm. 

I0  =  -2B(3)log  - -  +  DO) 


(87) 


•i-l 


where  D(p)  = 


i  2 

■rtf 


1-X‘ 


2 


log(m^-E-pcos0)  - ^  dx 


(l-3x) 


2 


-1 


(88) 


D(p)  is  a  well  behaved  quantity.  The  explicit  evaluation  of  it  is  not 
necessary,  because  x;hen  we  take  the  ratio  of  the  spectra  it  will  drop  out. 


n 


Tl 

iIBl 


2 


t  | 2 

°!  23  2W  26 


r  mF6E-m  m  . 

2  -2B(3)log  —  ----- 1 T-J  +  CO)  +  DO) 


(89) 


Now,  let  us  evaluate  the  cross  product  term  i.e.  eq.  (77). 


k-jS 


=  |U(~^)  (W)(1-y5)V  |v(l+iY5)^U 


:Z,Tr 

X 


1  /P.£ 


X 


2ymv  p 


("k*-)  (  ^  ^ 


=E-X-  (2^)  " 

X  P*k 


(p.k)(q,f  k)-(p.q)(K.£^)  +  (p.£^)(K.q)  + 


L 


. 


' 
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+(p.k)(q.£^')-(p.q)(k.£^)+(p.£A-)(k.q)  * 

=S-2—  (£47)h(p.K)(q.£A0-t2(p.£A-)(K.q)-(p.p)(q.£A') 

X  p.k  - 

-(P.k')(q.£i)-(p.£i)(q.k')]  (90) 


In  which  we  have  made  use  of  the  equality  K+k  =  2K-p-q-k!  and  the  con¬ 
ditions  Ko£^=0  and  k,e£^=0  ; 

and  have  neglected  the  term  which  is  multiplied  by  (q.q)  =  m^  since 
it  is  quadratic  in  mv,  and  small  compared  to  the  other  terms.  Substitut¬ 
ing  eq.  (90)  back  in  the  expression  for  ,  we  obtain 

ID 


*1*  = e 


2  :l2 


2 


(2ttP 


m  ,m 


P _ 

p.k' 


V  1 


(q.£A)'  2(p.K)-(p.p)-(p.kI  )  +  (p.6A) 


-  2(K.q)-(q.k') :  (2rr)  ^(K-k'-p-k-q) 


(91) 


*  ttfe  have  ornitted  the  terms  like  pU4iyq^£a0  Note  that 

c°  has  only  the  space  components, 

i.e.  1°  =  0,  f1  ^  0  i  =  1,2,3. 

ci 

Expanding  the  tensor  terms  one  will  get  terms  like  c  . . . c  _ 

In  any  cartesian  frame  one  will  have  b1  with  angles  0,  4  for  £. 
Therefore,  one  will  get  quantities  like 

)  ’cosG,  )  sinGcos'l',  and  y  sinGsin^. 

But  the  average  of  all  of  them  is  zero. 

tt  tt2tt  rr2rr 

/cosGdG  =0,  f  /sinGcosYdGdf  =0,  /  /sinGsintdGdT  =  0. 

o  00  00 

Hence,  these  terras  will  not  contribute. 


.  >-  • 


J 
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When  summing  and  averaging  over  the  polarization  vectors  of  the  photon, 
we  have 


(p.£ 


1 

2 


(p.6X)(q.£X)  =  i 


 (p.k* )(q.k» ) 


p.q  - 


And  we  also  have  the  identities: 


p.k'=  Ecj*  -  p.k’ 
p.K  =  m^E 


p.p  =  m. 


q.k1  =  £w*  -  q.k1 
q.K  = 


Putting  all  these  quantities  into  eq.  (91) »  we  get 


X  _  e2Q'2  1 

IB  =  (ar)3  V* 


( 2mgE-m^-Ecj'  +p . k 1 ) ( 


2  1 p.k1 | 2 


)  + 


+  (Em^-W+q.k'Xp.q-  (2n)4S(K-k'-p-k-q) 

p  #  K 1 


k'2 


(92) 


Let  the  angle  between  p  and  k’  be  0,  the  angle  between  q  and  k1  be  0’ 

and  the  angle  between  p  and  q  be  If  we  take  the  direction  of  p  as 

the  direction  of  the  z-axis  and  q  to  be  in  the  x-z  plane,  then  the  re¬ 
lationship  between  0  and  0’  is  given  by 


cos0’  =  cos0cos^  -i-  sin0cosl]Jsint/) 


(93) 


In  eq.  (93)  f  is  the  azimuthal  angle 


. 


U 4 


z 


Fig. 
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Recall  that  in  calculating  the  decay  probability,  we  have  to 


integrate  with  respect  to 


d3p  A  d\  d3k' 
2E  2oj  26  2w* 


For  the  moment,  we 


leave  out  the  integration  with  respect  to  the  first  three  variables. 


and  do  only  the  integration  with  respect  to 


d\» 

2u>» 


I Y  j  d3k» 

1  e2G'2  1  \ 

1 

-  — 

(  O™  T  T?  «  .itTli  7T  J  ! 

IB  2 co*  - 

\  (2tt)3  V»  (  (EW- 

P 

k» 

cosG) 

V 4CIHg-Xl>-»ni  +  p  1 1  K"  |  COStf,/ 

;  p,q»  j  p  1 1  q  !  (cos^Gcos'FfsinG  cosG  cosl)j  si nf)  '  + 

\ 

|  am  ,  -r-  1  «wsr  '  *3  '"> 

+i  2mr/(-=(.Ldf  + ! q i  1  k*  j  (cosGcos^+sinGcos^sin^ )  ipl  sin  0 

L  ^ 

llAkld^dcosO  (2rr)46(K»k«»p-k-q) 

2co’ 


>  1- , 
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e2G'2  1 


(2tt)4  2>mw 


Kax  |  (2mKE-m2  )  I  <ll  cosf-2mK«  I  P 


-  (6|i+iog  £g)!  -ik^jp.q 


L3  “  f  (‘23  +  1°g  1^} 


jp  (log  i±|  -  |g3-2) 


+  P 


£ 


log  -  -^(63+log  ^|) 
L  1-P  ^  l-P_ 


j 


(2tt)  6(K-k»-p-k-q) 


(94) 


We  can  see  from  eq.  (94)  that  the  correction  due  to  the  cross-term  is 
well  behaved  for  all  muon  momenta.  Let  us  write  the  cross-term  as: 


i*i  £l  £i  ii  £*l  =  A(@)£p£k£q 

IB  2E  2u  26  2o'  2E  2u  26 


(95) 


The  evaluation  of  the  integral 
the  two  previous  cases. 


y 

T^gj  is  much  simpler  than 


ipti 

IB 


2  e  G 


2r;i2 


(2tt)^  sX 


2,0  \4 


U^(l-iy^)vJ'c(2n)  6(K-k*-p-k-q) 


e2G>2  1  v-  T/  /  x^//,  .  x  ,o  x4. 
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where  A"  = 
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Finally,  the  radiative  correction  due  to  the  real  photon 
emission  process  is 
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Note  that  the  infrared  divergent  term  which  arises  here  is  equal  but 
opposite  in  sign  to  that  from  the  virtual  photon  process  corrections 
(c.f.  eq,  (68)  Chapter  III). 
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CHAPTER  V  THE  SPECTRA  OF  MUON  AND  PION 


5.1  THE  SPECTRUM  OF  THE  MUON 


When  evaluating  the  muon  spectrum  we  carry  out  all  integra¬ 
tions  except  the  one  over  the  muon  momentum  in  the  expression  for  the 
decay  probability*  Moreover,  the  virtual  photon  correction  and  the 
real  photon  correction  contribute  incoherently  to  the  spectrum,  because 
they  lead  to  final  states  with  different  number  of  particles. 

The  radiative  corrections  from  the  virtual  photon  processes 
to  the  muon  spectrum  are 


dP1 (E,m^)  = 
From  eq.  (68) 
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where  we  have  put  2(l-m^EFQ) = 2B(3).  Now,  we  take  the  square  of  the 
modulus  of  T  and  sum  over  the  spin  states  of  neutretto,  we  get 
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Where  TQ  is  the  uncorrected  decay  amplitude,  which  was  defined  earlier 
(c.f.  eq.  (67)  Chapter  III),  and  R  is  the  contribution  from  the  cross- 
product  terms. 


R  =  2Be2>'  ^  |_U(2K-m;J)(l-ir5)V+D(.mv)(l+ir5)V_ 


j_V(l+iY5)(f2+f,»t)U  +  V(l-ir5)l&iyDj  (2tt)  ‘'d(K-k-p-q)  (99) 

In  both  eq.  (98)  and  eq.  (99)*  terms  of  order  higher  than  e  have  been 
neglected.  Expanding  the  expression  for  R,  we  have 
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p  p 

The  rest  of  the  R's  are  of  order  m^  in  comparison  to  R.  Since  m^  is 
small,  we  may  neglect  them  and  keep  only  R-^. 

Hence,  the  muon  spectrum  with  virtual  photon  correction  will  be 
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The  real  photon  correction  to  the  muon  spectrum  is 
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It  follows  from  eq.  (97)»  that 
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The  final  spectrum  after  radiative  corrections  is 


dP(E,mu)  =  dP'CE.m^)  +  dP"(E,m„) 


50 


X-|(2BO)log(mK6E-miTml))-CO)-DO)+if1)]  +R  + 


|  2C  2to 


)  d^q  d^k 
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In  eq.  (103)  the  explicit  infrared  divergent  terms  have  cancelled  but 
there  remains  a  residual  divergent  term  log(m^6E-m^mv>)  which  tends  to  -<» 
as  the  end-point  of  the  spectrum  is  approached.  This  reflects  the  lack 
of  cancellation  of  infrared  divergences. 


The  muon  spectrum  is  now  obtained  by  carrying  out  the  q  and  k 
integrations.  One  of  these  integrations,  say  the  q-integration,  is 
trivial  because  of  the  presence  of  momentum  conserving  delta-function. 
The  range  of  the  last  integration  variable  is  limited  by  energy  momentum 
conservation.  We  need  not  know  this  range  for  the  calculation  of  the 
muon  spectrum  because  the  important  terms  in  eq.  (103)  are  all  k-inde- 
pendent  and  we  shall  only  be  interested  in  the  ratio  dP(E,m^)/dP(E,0) . 
however  for  the  sake  of  completeness  we  shall  digress  here  for  a  moment 
and  evaluate  the  range  of  energy  available  to  the  pion  when  the  muon 
energy  is  fixed.  The  calculation  of  the  range  of  muon  energy  for  a 
fixed  value  of  pion  energy  proceeds  in  an  exactly  analogous  fashion. 

In  the  K-meson-rest-frame , 
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or 


k 

max 

min 


(105) 


We  shall  carry  out  the  calculation  for  the  case  m^sO  only.  Then 
the  energy  of  the  neutretto.  By  energy  conservation 

6  =  mj£  -  E  -cj  (106) 

Dropping  the  indices  'max1  and  ’min'  and  substituting  eq.  (106)  in  eq. 
(105)  we  get 
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In  an  analogous  fashion  for  every  pion  energy  cj,  the  energy  range  of  the 
muon  is 
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We  can  now  show  from  eq.  (107)  that  if  the  muon  is  observed  near  its 
maximum  energy  then  the  pion  energy,  cj ,  has  a  very  limited  range  close 
to  its  maximum  value.  The  absolute  energy  maxima  are 
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2  2  2 
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Suppose  we  observe  muons  with  energy 
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max 
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and  momentum 

P  =  Praax  "  6P 

(113) 

with 

*P  -  5E/iW 

(114) 

Substituting  eq. 

(112)  -  eq.  (114)  in  eq.  (107)  we 

find  that  the  minimum 

pion  energy  is 

"  =  Haax  -  *  « 

Hmax 

(115) 

where  wmax  is  given  by  eq.  (110)  and  l/3max  >  1.  Thus  the  pion  carries 
away  almost  its  maximum  energy  when  the  muon  does  so  too. 


We  will  now  return  to  the  discussion  of  the  muon  spectrum. 
The  R  term  in  eq.  (103)  contains  terms  fg  and  f^.  These  terms  as  men¬ 
tioned  in  Chapter  III,  are  small  compared  with  the  terms  which  diverge 
at  the  end-point  of  the  muon  spectrum.  The  A  terms  as  we  can  see  from 
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eq.  (9*0  and  eq.  (96)  of  Chapter  IV,  are  well  behaved  and  do  not  give 
rise  to  any  divergent  term.  We  may  drop  all  these  terms  in  the  expres¬ 
sion  of  the  muon-spectrum  near  the  end-point.  Therefore  the  muon- 
spectrum  near  the  end-point  is 
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The  ratio  dP(Em,my)/dP(Em,0)  can  now  be  written  as  follows: 
dP(Em,mi;)/dP(Em,0) 
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where  dPo(Em,m^)  and  dPQ(Em,0)  are  the  uncorrected  muon  spectrum  with 
m^O  and  that  with  m^=0  respectively.  Their  expressions  were  given  in 
eq.  (6)  and  eq.  (7)  of  Chapter  I.  The  value  of  (3^  with  m^O  differs 
from  with  rriy=0  by  terms  of  order  of  m^,.  We  may  take  the  functions 
C(j3  ),  b(j3m)  and  f-^  to  be  equal  in  both  cases.  Then,  we  can  write  eq. 
(117)  as 
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The  ratio  of  the  uncorrected  spectra  dP  (E  .m^/dP^E  ,0)  is  the  same 
as  that  in  eq.  (12)  of  Chapter  I.  Hence,  eq.  (118)  becomes 
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5.2  THE  SPECTRUM  OF  THE  PION 


In  evaluating  the  full  pion  spectrum  (i.e.  the  spectrum  for 


all  pion  energies)  we  run  into  difficulties.  These  difficulties  have 
their  origin  in  the  fact  that  the  radiative  correction  terms  are  func¬ 
tions  of  the  muon  momentum  and  an  integration  over  the  muon  spectrum 
has  to  be  done  to  get  the  pion  spectrum.  Though  a  numerical  integra¬ 
tion  is  feasible,  it  is  extremely  difficult,  if  at  all  possible,  to 
get  an  analytic  form  for  the  spectrum.  One  can  get  around  this  dif- 
iculty  if  one  is  interested  in  only  the  end-point  (and  its  neighbour¬ 
hood)  of  the  pion  spectrum.  In  this  case  the  muon  energy  is  close  to 
its  maximum  as  was  shown  in  the  discussion  on  energy  ranges  in  the 
section  5.1  on  the  muon  spectrum.  The  p-dependent  integrand  can  be 

treated  as  a  constant  with  p  replaced  by  p  . 

max. 


The  calculation  proceeds  in  an  analogous  fashion  to  that  of 


the  muon  spectrum.  We  only  have  to  express  the  maximum  soft  photon 
energy  in  terms  of  the  pion  energy  rather  than  the  muon  energy.  The 
precise  relation  is 
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The  pion  spectrum  will  then  be  (analogous  to  eq.  (116)) 
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±^y  following  the  same  reasoning  as  that  for  the  muon  spectrum  we  obtain 
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5.3  ENERGY  CORRELATION  BETWEEN  MUON  AND  PION 


Let  us  consider  the  phase  space  integral  with  radiative  cor¬ 


rections  , 


To|2[1*S[2B(^log(vE-v,i')-c(@)-D(@)+*fi}] (123) 


In  eq.  (123)  we  have  dropped  the  term  R,  A  and  A"  because  we  are  going 
to  look  at  the  correlation  at  the  high  energy  end,  where  the  contribu¬ 
tion  from  these  terms  is  insignificant.  Eq.  (123)  can  be  written  as 
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where  >7  is  the  angle  between  pion  and  muon  momenta.  The  integration 

/> 

over  the  angle  6 1  cos  ^  -( ♦ . .  )j  dcosr}  is  equal  to  1.  Eq.  (124)  becomes 
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where  p(E,cj)  is  the  energy  correlation  between  the  pion  and  muon  for  the 
particular  model  namiltonian  we  have  chosen. 
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The  last  2  terms  in  the  eq.  (126)  when  compared  with  the  first  term  are 
2  4 

of  order  my  and  my  respectively  and  so  are  negligible. 
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Since  (p+q)  =  (K-k)  ,  if  the  K  meson  decays  at  rest,  then 
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In  deriving  eq.  (128)  we  have  made  use  of  eq.  (110)  and  have  neglected 
terms  of  order  m^.  Substitute  eq.  (128)  into  eq.  (127),  we  get 

T0|2  =  G'2m^J  |  8L2E£  ■  mK(wmax  "  ">]  "  8mKm?£  +  2n^mK(<Jma3c-  u>)J  (129) 


V(  i  v  f  <- j,  :,J|  >  ,  !  ;*<  *  •  ,CYX-X>vj 


■ 


**'■ 


59 


Finally,  the  energy  correlation  between  pion  and  muon  is 
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Taking  the  ratio  of  p(E,(J)  with  m^O  to  that  with  m^=0,  we  have 
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We  can  see  from  eq.  (129)  that 


is  a  function  of  E  and  6w=6J»  -w. 
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The  only  place  in  this  ratio  where  mv  comes  in  is  through  the  cosr? 
integration  and  the  Innerbrerasstrahlung  correction,  i.e.  logCm^SE-nyn^) . 
The  value  of  cost]  in  the  case  m^O  differs  from  the  value  in  the  case 
niysO  by  terms  quadratic  in  m^  and  consequently  we  do  not  expect  a  large 
effect  from  the  cosh  term.  The  only  important  term  appears  to  be  the 
innerbremsstrahlung  correction  term  log(mK6E-m  ra^).  Recall  that  when 
E=Em-6E,  GJ  is  very  close  to  (Jm  and  we  can  now  look  at  p(E,o>)  as  a  func¬ 
tion  of  E  for  a  particular  value  of  cj,  say  Hence,  the  ratio  of  the 

energy  correlations,  for  a  particular  value  of  pion  energy,  near  the 
high  muon  energy  end  is 
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The  reason  for  dropping  the  terms  C(bm),  D(0m)  and  f^»  and  for  taking 
12 


the  two 


to  be  equal  is  the  same  as  that  for  the  muon  spectrum 


case  (Section  5.1). 


i  o: 
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CHAPTER  VI  CONCLUSION  AND  DISCUSSION 


The  results  of  the  calculations  performed  in  the  present  work 
are  displayed  in  eq.  (119) »  ©q.  (122)  and  eq.  (132)  of  Chapter  V.  We 
see  from  these  equations  that  the  radiative  corrections  to  the  ratio  of 
the  spectra  with  and  without  the  neutretto  mass  comes  in  the  form  of  the 
logarithmic  term 
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or 
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|  BO  )log(l~t!_l). 
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This  term  was  referred  to  as  the  universal  logarithmic  singularity  in 

(3) 

the  paper  by  Allcock  (1965) ,  The  presence  of  this  terra  reflects  the 
lack  of  cancellation  of  the  infrared  divergence.  If  we  observe  the 
spectrum  of  one  of  the  decay  particles,  say  the  spectrum  of  the  muon, 
near  the  high  energy  end,  we  expect  the  ratio  dP(E,my)/dP(E,0)  to  be 
deviated  from  that  without  corrections  by  an  amount 
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When  6E  — —  ,  the  logarithmic  term  tends  to  -  «>.  But  a  measurement 
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of  the  muon  spectrum  in  the  region  where  6E  is  very  close  to  m^/m^m^ 
would  yield  an  upper  limit  of  neutretto  mass.  The  upper  limit  of  my 
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can  be  set  by  looking  at  the  deviation  from  unity  of  the  graph  of  the 
ratio  dP(E,m^)/dP(E,0) ,  plotted  against  the  energy  of  the  muon,  near 
the  end-point*  The  dP(E,my)  will  be  given  directly  from  experiment, 
while  dP(E,0)  can  be  calculated  theoretically.  Without  the  radiative 
corrections  the  ratio  would  be  essentially  the  phase  space  ratio 
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and  that  has  been  fully  discussed  in  both  the  papers  by  Dennery  and 
Primakoff  (1963)^^  and  by  Ginsberg  (1965)^^ .  But  if  B(0m)  is  not  too 
small,  the  contribution  from  the  radiative  corrections  may  become  sig¬ 
nificant.  In  the  KIJO  decay,  mv  »  m  or  mu,  0  1,  thus  B(0  )  will  be 

rj  iv  rr  r  m  m 

of  the  order  of  unity,  so  we  expect  the  radiative  corrections  to  be 
quite  significant.  For  the  muon  spectrum,  the  correction  to  the  ratio 
is  1$  for  6E  a*  1.2m  m^/m^  or  0.33  My  In  the  case  of  the  pion  spectrum, 
because  of  the  presence  of  the  term  m^,my  as  against  the  term  rn_my  in  the 
muon  spectrum  case  we  shall  have  to  get  closer  to  the  end-point  before 
we  can  see  a  significant  departure  from  unity  in  the  ratio.  For  in¬ 
stance,  the  correction  is  1$  for  dco  &  1.2m^m^/m^.  or  0.25  m^.  Moreover 
the  electrical  neutrality  of  pion  makes  its  spectrum  experimentally 
less  tractable  than  that  of  the  muon.  For  the  ratio  of  the  energy  cor¬ 
relations,  the  correction  will  be  essentially  the  same  as  that  in  the 
muon  spectrum  case.  However,  the  correlation  experiments  are  more  dif¬ 
ficult  to  perform  than  the  direct  spectrum  measurement. 
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APPENDIX  A  PHASE-SPACE  INTEGRATION  IN 
THE  TIME-LIKE-G-FRAMB 


We  have  in  Chapter  I  eq.  (2)  the  phase-space  integral  for  the 
muon  spectrum. 


dP(E,my)  = 


6(  G-k-q) 

2  U)  2£ 


(A.l) 


In  which  we  have  set  G=K-p.  Now,  let  us  write  — as: 

2£ 


^2.  =  d^qdq°6(q°  -q2-m2)  =  d^qdq°6(q°  -£2) 


(A.2) 


2  2  2 

where  €  =q  -Hti  .  Substituting  eq.  (A.2)  into  eq.  (A.l),  then  we  get 


r\ 


dP(E.m^)  = 


d4q  6(qo2-€2)64(G-k-q) 

£(j) 


(A.3) 


In  eq.  (A.3)  the  q  integration  is  trivial 


dP(E,mu) 


A 

26J 


(G-k)1 


(A.4) 


This  expression  is  Lorentz  invariant  because  both 


di  „  d\6 

2<J 


and 


(G-k) 


.  • .  y 


■  - 


ii 


are  Lorentz  invariant. 


So  we  can  work  in  the  time-like  G  frame  i.e.  G^  = ( 


where  eq.  (A. 4)  becomes 


dP(E,m.,)  = 


d^k 


+  k2  -  2G°k0  -  m2)j 


d^k  A /  r2  o r>Oi  ,  c  c\ 

~  o(  u  -  4'a  k  +  m  —  m, ,) 


.2  JZ' 


2. 


to 


IT 


UJ 


&  JL  6(k 

2(Jo  o 


p2  2 

Lr  +  m 
O  TT 


2G 


2 


-  -  oo 

we  can  write  kdk  =  k  dk  ,  so 


4-^  =  Ztt  j  kdk° 

Z(jJ 


Substituting  eq.  (A. 6)  into  eq.  (A.5)»  we  have 


dP(E,mJ  = 


2rr  Mil  6(k 
2G°  0 
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2  2  2 

Vy”1^ 

2G° 


=.-  rr 


evaluated  at  k  = 
o 


P2  2  2 
G  +m  -m.. 
o  tt  V 

2G° 


r2  ,2  2 

k  =  k  -  m 

O  TT 


fr2  2  212 

G  +m  -m  ! 


o  n  \>  i  2 

- :  -  m 

o  i  n 


2G 


0,0,0). 


(A,5) 


(A.6) 


(A. 7) 


... 
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_  1 


r,4  or2,  2  (  2  2x2 

Go-2Go(mrr  +D1^  +  ^Vmv) 


(A.8) 


Put  this  result  into  eq.  (A. 7)  to  obtain 


dP(S,mJ  s 


IT 

2  o  L 


r4  op2/  2  2n  ^  ,2  2x2 

-  2a0(mn+niy )  +  (m^-m^ ; 


2  2  2  2 
This  is  in  the  time-like  G  frame,  GQ  =  G  ,  so  one  can  replace  Gq  by  G 

all  along. 


dP(E,m^)  =  -% 
2G^ 


G4  -  2G2(m2+m2)  +  (m2-n£)2 
TT  v  TT  v 


(A.9) 
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APPENDIX  B 


INTEGRATION  OF 


2d\' 

2iA* 


IN 


SPHERICAL  COORDINATES 


To  evaluate  this  integral  we  assign  to  the  photon  a  small 
mass  X.  so  as  to  avert  the  divergence.  Then  co*  =k,fc  +  Ar  , 

Let  0  be  the  angle  between  p  and  k* 
p/E=  B  =  velocity  of  the  muon 

k!/ojs  =  U,  corresponds  to  the  velocity  of  the  photon. 


Upon  summing 
have 


p.k5 


over  the  polarization  vectors  of  the  photon,  we 


p.k5 


d3k5  i  k2  (1-152cos20)  d^k5 
!  (l-j3tfcos8)2 


(B.l) 


We  take  the  direction  of  the  p  as  the  direction  of  the  z  axis,  and  do 
the  integration  in  spherical  coordinates. 
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p.k* 


2  3 
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pi 


=  2rr 


sin0d0 


K5 

max 

P' 
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(l-lj2cos29)  2dk* 


(l-0UcOS0r  w 

We  can  replace  the  t-integration  by  the  V  integration  in  eq.  (B.2) 


(B.2) 
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Since  U  =  k' /(k'2+A.2)2  then  d  V  =  — ~- 


(k,2+X2p/2 


dk' 


0J'  1-tf2 


and  lf=0  when  k’=0,  1)=1-|  W' 


max 


when  k'=k' 

max 


Substituting  these  quantities  into  eq.(B.2)  and  letting  x=:cos0  we  obtain 
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+1  i-Hax 


dx 
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2  2 

Tixll  dlS 

(i  -mf 


=  I»  +  I" 


(B.3) 


The  integration  of  I’  and  I"  are  quite  tedious,  their  results  are  given 
as  follows: 


+1 
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vi 


where  the  L-f unctions  are  the  Spence  functions.  Finally  we  have  the 
expression 
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2d\' 

2uS« 


+1 


b-t  ? 

max  l_x 
K  (1-fix)2 


dx  -  4B(d)(log2-2)  + 


+  |[lO)-L(-@)]  -  i  [L(^2)-L(i=£) 


+  2  - 


-  2.  log  i=£  (log  +  2) 
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(B.6) 


